(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 10.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     38309,       1081]
NotebookOptionsPosition[     36331,       1010]
NotebookOutlinePosition[     36685,       1026]
CellTagsIndexPosition[     36642,       1023]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{"ClearAll", "[", "\"\<Global`*\>\"", "]"}]], "Input"],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
   "Can", " ", "be", " ", "used", " ", "for", " ", "Dirichlet", " ", "and", 
    " ", "traction", " ", "BCs"}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{"Desired", " ", "solution"}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"u", "[", 
      RowBox[{"x_", ",", "y_", ",", "t_"}], "]"}], ":=", 
     RowBox[{"2", "Pi", " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"2", " ", "Pi", " ", "x"}], "]"}], " ", 
      RowBox[{"Cos", "[", 
       RowBox[{
        RowBox[{"2", " ", "Pi", " ", "y"}], "-", 
        RowBox[{"2", " ", "Pi", " ", "t"}]}], "]"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"v", "[", 
      RowBox[{"x_", ",", "y_", ",", "t_"}], "]"}], ":=", " ", 
     RowBox[{
      RowBox[{"2", " ", "Pi", " ", 
       RowBox[{"Sin", "[", 
        RowBox[{"2", " ", "Pi", " ", "x"}], "]"}], 
       RowBox[{"Sin", "[", 
        RowBox[{
         RowBox[{"2", "Pi", " ", "y"}], "-", 
         RowBox[{"2", " ", "Pi", " ", "t"}]}], "]"}]}], " ", "+", " ", 
      RowBox[{"Sin", "[", 
       RowBox[{
        RowBox[{"2", " ", "Pi", " ", "x"}], "-", 
        RowBox[{"2", "Pi", " ", "t"}]}], "]"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"p", "[", 
      RowBox[{"x_", ",", "y_", ",", "t_"}], "]"}], ":=", 
     RowBox[{"2", "Pi", " ", 
      RowBox[{"Sin", "[", 
       RowBox[{
        RowBox[{"2", "Pi", " ", "x"}], "-", 
        RowBox[{"2", "Pi", " ", "t"}]}], "]"}], 
      RowBox[{"Cos", "[", 
       RowBox[{
        RowBox[{"2", "Pi", " ", "y"}], "-", 
        RowBox[{"2", "Pi", " ", "t"}]}], "]"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"d", "[", 
      RowBox[{"x_", ",", "y_"}], "]"}], ":=", 
     RowBox[{"0.1", "-", 
      RowBox[{"Sqrt", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{"x", "-", "0.5"}], ")"}], "^", "2"}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"y", "-", "0.5"}], ")"}], "^", "2"}]}], "]"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"\[Rho]", "[", 
      RowBox[{"x_", ",", "y_", ",", "t_"}], "]"}], ":=", 
     RowBox[{"\[Rho]0", "+", 
      RowBox[{
       RowBox[{"\[Rho]1", "/", "2"}], 
       RowBox[{"(", 
        RowBox[{"1", "+", 
         RowBox[{"Tanh", "[", 
          RowBox[{
           RowBox[{"d", "[", 
            RowBox[{"x", ",", "y"}], "]"}], "/", "\[Delta]"}], "]"}]}], 
        ")"}]}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"\[Mu]", "[", 
      RowBox[{"x_", ",", "y_", ",", "t_"}], "]"}], ":=", 
     RowBox[{
      RowBox[{"\[Mu]1", " ", 
       RowBox[{"Sin", "[", 
        RowBox[{"2", " ", "Pi", " ", "x"}], "]"}], 
       RowBox[{"Cos", "[", 
        RowBox[{"2", " ", "Pi", " ", "y"}], "]"}]}], "+", "\[Mu]1", " ", "+", 
      "\[Mu]0"}]}], ";"}], "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.713710130967271*^9, 3.71371016525543*^9}, {
   3.713710201737459*^9, 3.713710227808384*^9}, {3.713711199666253*^9, 
   3.7137112047861977`*^9}, {3.713714563101632*^9, 3.7137146175729103`*^9}, {
   3.7137148887951717`*^9, 3.713714941641158*^9}, {3.713715050278842*^9, 
   3.713715081235734*^9}, 3.7137151226289988`*^9, {3.713715193003085*^9, 
   3.713715194252469*^9}, {3.713719829770111*^9, 3.713719885616111*^9}, {
   3.716904649871188*^9, 3.7169047034396467`*^9}, {3.716904756071067*^9, 
   3.716904763841667*^9}, {3.7173540928966227`*^9, 3.717354122818575*^9}, {
   3.717790754816205*^9, 3.71779077247313*^9}, {3.717791239840149*^9, 
   3.717791239986281*^9}, {3.7177912715201893`*^9, 3.7177912716002398`*^9}, {
   3.7177913104170847`*^9, 3.7177913630227003`*^9}, {3.717806555069543*^9, 
   3.717806558832151*^9}, {3.7295151681376534`*^9, 3.729515169263995*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"tau11", "[", 
    RowBox[{"x_", ",", "y_", ",", "t_"}], "]"}], "=", "  ", 
   RowBox[{"D", "[", 
    RowBox[{
     RowBox[{"u", "[", 
      RowBox[{"x", ",", "y", ",", "t"}], "]"}], ",", "x"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"tau12", "[", 
    RowBox[{"x_", ",", "y_", ",", "t_"}], "]"}], " ", "=", " ", 
   RowBox[{
    RowBox[{"1", "/", "2"}], 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"D", "[", 
       RowBox[{
        RowBox[{"u", "[", 
         RowBox[{"x", ",", "y", ",", "t"}], "]"}], ",", "y"}], "]"}], "+", 
      " ", 
      RowBox[{"D", "[", 
       RowBox[{
        RowBox[{"v", "[", 
         RowBox[{"x", ",", "y", ",", "t"}], "]"}], ",", "x"}], "]"}]}], 
     ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"tau22", "[", 
    RowBox[{"x_", ",", "y_", ",", "t_"}], "]"}], "=", "  ", 
   RowBox[{"D", "[", 
    RowBox[{
     RowBox[{"v", "[", 
      RowBox[{"x", ",", "y", ",", "t"}], "]"}], ",", "y"}], "]"}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.7173541393691998`*^9, 3.71735415067225*^9}, {
  3.717791265564101*^9, 3.717791266924889*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"Divergence", " ", 
    RowBox[{"free", "?"}]}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"D", "[", 
    RowBox[{
     RowBox[{"u", "[", 
      RowBox[{"x", ",", "y", ",", "t"}], "]"}], ",", "x"}], "]"}], "+", 
   RowBox[{"D", "[", 
    RowBox[{
     RowBox[{"v", "[", 
      RowBox[{"x", ",", "y", ",", "t"}], "]"}], ",", "y"}], 
    "]"}]}]}]], "Input",
 CellChangeTimes->{{3.713714627201284*^9, 3.713714643136862*^9}, {
  3.7137150370093937`*^9, 3.713715041758703*^9}, {3.7137150838464527`*^9, 
  3.713715084650531*^9}, {3.717791345063168*^9, 3.717791345705038*^9}}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{
  3.713714644055319*^9, 3.713714948053357*^9, 3.713715085121151*^9, 
   3.713715144571857*^9, 3.7137151962065907`*^9, 3.71371544182511*^9, 
   3.713719900407336*^9, 3.713720451520116*^9, 3.7156965043275337`*^9, 
   3.7159696431675034`*^9, 3.7169047077321577`*^9, 3.7169047658532143`*^9, 
   3.7173541616910067`*^9, 3.717358317990509*^9, 3.717448734480977*^9, 
   3.717791081084073*^9, 3.717791268619244*^9, {3.717791330038245*^9, 
   3.71779136640241*^9}, 3.717791615519669*^9, 3.71780656220935*^9, 
   3.717857650839448*^9, 3.7185455346465282`*^9, 3.729466252608985*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"Required", " ", "forcing", " ", "function"}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"fx", "[", 
     RowBox[{"x_", ",", "y_", ",", "t_"}], "]"}], "=", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"\[Rho]", "[", 
       RowBox[{"x", ",", "y", ",", "t"}], "]"}], 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{"u", "[", 
           RowBox[{"x", ",", "y", ",", "t"}], "]"}], ",", "t"}], "]"}], " ", 
        "+", " ", 
        RowBox[{
         RowBox[{"u", "[", 
          RowBox[{"x", ",", "y", ",", "t"}], "]"}], 
         RowBox[{"D", "[", 
          RowBox[{
           RowBox[{"u", "[", 
            RowBox[{"x", ",", "y", ",", "t"}], "]"}], ",", "x"}], "]"}]}], 
        " ", "+", " ", 
        RowBox[{
         RowBox[{"v", "[", 
          RowBox[{"x", ",", "y", ",", "t"}], "]"}], " ", 
         RowBox[{"D", "[", 
          RowBox[{
           RowBox[{"u", "[", 
            RowBox[{"x", ",", "y", ",", "t"}], "]"}], ",", "y"}], "]"}]}]}], 
       ")"}]}], " ", "+", " ", 
     RowBox[{"D", "[", 
      RowBox[{
       RowBox[{"p", "[", 
        RowBox[{"x", ",", "y", ",", "t"}], "]"}], ",", "x"}], "]"}], " ", "-",
      " ", 
     RowBox[{"2", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"\[Mu]", "[", 
            RowBox[{"x", ",", "y", ",", "t"}], "]"}], " ", 
           RowBox[{"tau11", "[", 
            RowBox[{"x", ",", "y", ",", "t"}], "]"}]}], ",", "x"}], "]"}], 
        "+", 
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"\[Mu]", "[", 
            RowBox[{"x", ",", "y", ",", "t"}], "]"}], " ", 
           RowBox[{"tau12", "[", 
            RowBox[{"x", ",", "y", ",", "t"}], "]"}]}], ",", "y"}], "]"}]}], 
       ")"}]}]}]}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"fy", "[", 
     RowBox[{"x_", ",", "y_", ",", "t_"}], "]"}], "=", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"\[Rho]", "[", 
       RowBox[{"x", ",", "y", ",", "t"}], "]"}], 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{"v", "[", 
           RowBox[{"x", ",", "y", ",", "t"}], "]"}], ",", "t"}], "]"}], " ", 
        "+", " ", 
        RowBox[{
         RowBox[{"u", "[", 
          RowBox[{"x", ",", "y", ",", "t"}], "]"}], 
         RowBox[{"D", "[", 
          RowBox[{
           RowBox[{"v", "[", 
            RowBox[{"x", ",", "y", ",", "t"}], "]"}], ",", "x"}], "]"}]}], 
        " ", "+", " ", 
        RowBox[{
         RowBox[{"v", "[", 
          RowBox[{"x", ",", "y", ",", "t"}], "]"}], " ", 
         RowBox[{"D", "[", 
          RowBox[{
           RowBox[{"v", "[", 
            RowBox[{"x", ",", "y", ",", "t"}], "]"}], ",", "y"}], "]"}]}]}], 
       ")"}]}], " ", "+", " ", 
     RowBox[{"D", "[", 
      RowBox[{
       RowBox[{"p", "[", 
        RowBox[{"x", ",", "y", ",", "t"}], "]"}], ",", "y"}], "]"}], " ", "-",
      " ", 
     RowBox[{"2", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"\[Mu]", "[", 
            RowBox[{"x", ",", "y", ",", "t"}], "]"}], " ", 
           RowBox[{"tau12", "[", 
            RowBox[{"x", ",", "y", ",", "t"}], "]"}]}], ",", "x"}], "]"}], 
        "+", 
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"\[Mu]", "[", 
            RowBox[{"x", ",", "y", ",", "t"}], "]"}], 
           RowBox[{"tau22", "[", 
            RowBox[{"x", ",", "y", ",", "t"}], "]"}]}], ",", "y"}], "]"}]}], 
       ")"}]}]}]}]}]}]], "Input",
 CellChangeTimes->{{3.713715089011114*^9, 3.7137151124079742`*^9}, {
  3.713715166949651*^9, 3.7137151831024303`*^9}, {3.713715295717227*^9, 
  3.7137154214025183`*^9}, {3.713718298670486*^9, 3.7137183029435883`*^9}, {
  3.7137190607757874`*^9, 3.713719082418027*^9}, {3.71371912144843*^9, 
  3.713719136736127*^9}, {3.717354182465254*^9, 3.717354271509441*^9}, {
  3.7173543325691147`*^9, 3.717354346587216*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"4", " ", 
   SuperscriptBox["\[Pi]", "2"], " ", 
   RowBox[{"Cos", "[", 
    RowBox[{
     RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
     RowBox[{"2", " ", "\[Pi]", " ", "x"}]}], "]"}], " ", 
   RowBox[{"Cos", "[", 
    RowBox[{
     RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
     RowBox[{"2", " ", "\[Pi]", " ", "y"}]}], "]"}]}], "-", 
  RowBox[{"2", " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "8"}], " ", 
      SuperscriptBox["\[Pi]", "3"], " ", "\[Mu]1", " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"2", " ", "\[Pi]", " ", "x"}], "]"}], " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"2", " ", "\[Pi]", " ", "y"}], "]"}], " ", 
      RowBox[{"Cos", "[", 
       RowBox[{
        RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
        RowBox[{"2", " ", "\[Pi]", " ", "y"}]}], "]"}], " ", 
      RowBox[{"Sin", "[", 
       RowBox[{"2", " ", "\[Pi]", " ", "x"}], "]"}]}], "-", 
     RowBox[{"8", " ", 
      SuperscriptBox["\[Pi]", "3"], " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"2", " ", "\[Pi]", " ", "x"}], "]"}], " ", 
      RowBox[{"Cos", "[", 
       RowBox[{
        RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
        RowBox[{"2", " ", "\[Pi]", " ", "y"}]}], "]"}], " ", 
      RowBox[{"(", 
       RowBox[{"\[Mu]0", "+", "\[Mu]1", "+", 
        RowBox[{"\[Mu]1", " ", 
         RowBox[{"Cos", "[", 
          RowBox[{"2", " ", "\[Pi]", " ", "y"}], "]"}], " ", 
         RowBox[{"Sin", "[", 
          RowBox[{"2", " ", "\[Pi]", " ", "x"}], "]"}]}]}], ")"}]}], "-", 
     RowBox[{"2", " ", 
      SuperscriptBox["\[Pi]", "2"], " ", "\[Mu]1", " ", 
      RowBox[{"Cos", "[", 
       RowBox[{
        RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
        RowBox[{"2", " ", "\[Pi]", " ", "x"}]}], "]"}], " ", 
      RowBox[{"Sin", "[", 
       RowBox[{"2", " ", "\[Pi]", " ", "x"}], "]"}], " ", 
      RowBox[{"Sin", "[", 
       RowBox[{"2", " ", "\[Pi]", " ", "y"}], "]"}]}]}], ")"}]}], "+", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "8"}], " ", 
      SuperscriptBox["\[Pi]", "3"], " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"2", " ", "\[Pi]", " ", "x"}], "]"}], " ", 
      SuperscriptBox[
       RowBox[{"Cos", "[", 
        RowBox[{
         RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
         RowBox[{"2", " ", "\[Pi]", " ", "y"}]}], "]"}], "2"], " ", 
      RowBox[{"Sin", "[", 
       RowBox[{"2", " ", "\[Pi]", " ", "x"}], "]"}]}], "-", 
     RowBox[{"4", " ", 
      SuperscriptBox["\[Pi]", "2"], " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"2", " ", "\[Pi]", " ", "x"}], "]"}], " ", 
      RowBox[{"Sin", "[", 
       RowBox[{
        RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
        RowBox[{"2", " ", "\[Pi]", " ", "y"}]}], "]"}]}], "+", 
     RowBox[{"4", " ", 
      SuperscriptBox["\[Pi]", "2"], " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"2", " ", "\[Pi]", " ", "x"}], "]"}], " ", 
      RowBox[{"Sin", "[", 
       RowBox[{
        RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
        RowBox[{"2", " ", "\[Pi]", " ", "y"}]}], "]"}], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", 
         RowBox[{"Sin", "[", 
          RowBox[{
           RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
           RowBox[{"2", " ", "\[Pi]", " ", "x"}]}], "]"}]}], "-", 
        RowBox[{"2", " ", "\[Pi]", " ", 
         RowBox[{"Sin", "[", 
          RowBox[{"2", " ", "\[Pi]", " ", "x"}], "]"}], " ", 
         RowBox[{"Sin", "[", 
          RowBox[{
           RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
           RowBox[{"2", " ", "\[Pi]", " ", "y"}]}], "]"}]}]}], ")"}]}]}], 
    ")"}], " ", 
   RowBox[{"(", 
    RowBox[{"\[Rho]0", "+", 
     RowBox[{
      FractionBox["1", "2"], " ", "\[Rho]1", " ", 
      RowBox[{"(", 
       RowBox[{"1", "+", 
        RowBox[{"Tanh", "[", 
         FractionBox[
          RowBox[{"0.1`", "\[VeryThinSpace]", "-", 
           SqrtBox[
            RowBox[{
             SuperscriptBox[
              RowBox[{"(", 
               RowBox[{
                RowBox[{"-", "0.5`"}], "+", "x"}], ")"}], "2"], "+", 
             SuperscriptBox[
              RowBox[{"(", 
               RowBox[{
                RowBox[{"-", "0.5`"}], "+", "y"}], ")"}], "2"]}]]}], 
          "\[Delta]"], "]"}]}], ")"}]}]}], ")"}]}]}]], "Output",
 CellChangeTimes->{{3.7137151835747633`*^9, 3.713715197281507*^9}, 
   3.713715377068603*^9, {3.713715421668744*^9, 3.713715442688013*^9}, {
   3.713718300209652*^9, 3.713718303371324*^9}, {3.713719094128516*^9, 
   3.7137191371153383`*^9}, 3.7137204530821238`*^9, 3.715696505354093*^9, 
   3.715969644691567*^9, 3.716904709351152*^9, 3.716904766902685*^9, 
   3.717354248157126*^9, {3.71735432615724*^9, 3.717354346852295*^9}, 
   3.717358321406577*^9, 3.71744873648444*^9, 3.717791082128022*^9, 
   3.71779136874577*^9, 3.7177916176324883`*^9, 3.7178065650079308`*^9, 
   3.717857654516314*^9, 3.7294662537271423`*^9, 3.72946667576934*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"-", "4"}], " ", 
   SuperscriptBox["\[Pi]", "2"], " ", 
   RowBox[{"Sin", "[", 
    RowBox[{
     RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
     RowBox[{"2", " ", "\[Pi]", " ", "x"}]}], "]"}], " ", 
   RowBox[{"Sin", "[", 
    RowBox[{
     RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
     RowBox[{"2", " ", "\[Pi]", " ", "y"}]}], "]"}]}], "-", 
  RowBox[{"2", " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"2", " ", 
      SuperscriptBox["\[Pi]", "2"], " ", "\[Mu]1", " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"2", " ", "\[Pi]", " ", "x"}], "]"}], " ", 
      RowBox[{"Cos", "[", 
       RowBox[{
        RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
        RowBox[{"2", " ", "\[Pi]", " ", "x"}]}], "]"}], " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"2", " ", "\[Pi]", " ", "y"}], "]"}]}], "+", 
     RowBox[{"2", " ", 
      SuperscriptBox["\[Pi]", "2"], " ", 
      RowBox[{"(", 
       RowBox[{"\[Mu]0", "+", "\[Mu]1", "+", 
        RowBox[{"\[Mu]1", " ", 
         RowBox[{"Cos", "[", 
          RowBox[{"2", " ", "\[Pi]", " ", "y"}], "]"}], " ", 
         RowBox[{"Sin", "[", 
          RowBox[{"2", " ", "\[Pi]", " ", "x"}], "]"}]}]}], ")"}], " ", 
      RowBox[{"Sin", "[", 
       RowBox[{
        RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
        RowBox[{"2", " ", "\[Pi]", " ", "x"}]}], "]"}]}], "-", 
     RowBox[{"8", " ", 
      SuperscriptBox["\[Pi]", "3"], " ", "\[Mu]1", " ", 
      RowBox[{"Cos", "[", 
       RowBox[{
        RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
        RowBox[{"2", " ", "\[Pi]", " ", "y"}]}], "]"}], " ", 
      SuperscriptBox[
       RowBox[{"Sin", "[", 
        RowBox[{"2", " ", "\[Pi]", " ", "x"}], "]"}], "2"], " ", 
      RowBox[{"Sin", "[", 
       RowBox[{"2", " ", "\[Pi]", " ", "y"}], "]"}]}], "+", 
     RowBox[{"8", " ", 
      SuperscriptBox["\[Pi]", "3"], " ", 
      RowBox[{"Sin", "[", 
       RowBox[{"2", " ", "\[Pi]", " ", "x"}], "]"}], " ", 
      RowBox[{"(", 
       RowBox[{"\[Mu]0", "+", "\[Mu]1", "+", 
        RowBox[{"\[Mu]1", " ", 
         RowBox[{"Cos", "[", 
          RowBox[{"2", " ", "\[Pi]", " ", "y"}], "]"}], " ", 
         RowBox[{"Sin", "[", 
          RowBox[{"2", " ", "\[Pi]", " ", "x"}], "]"}]}]}], ")"}], " ", 
      RowBox[{"Sin", "[", 
       RowBox[{
        RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
        RowBox[{"2", " ", "\[Pi]", " ", "y"}]}], "]"}]}]}], ")"}]}], "+", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "2"}], " ", "\[Pi]", " ", 
      RowBox[{"Cos", "[", 
       RowBox[{
        RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
        RowBox[{"2", " ", "\[Pi]", " ", "x"}]}], "]"}]}], "-", 
     RowBox[{"4", " ", 
      SuperscriptBox["\[Pi]", "2"], " ", 
      RowBox[{"Cos", "[", 
       RowBox[{
        RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
        RowBox[{"2", " ", "\[Pi]", " ", "y"}]}], "]"}], " ", 
      RowBox[{"Sin", "[", 
       RowBox[{"2", " ", "\[Pi]", " ", "x"}], "]"}]}], "+", 
     RowBox[{"2", " ", "\[Pi]", " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"2", " ", "\[Pi]", " ", "x"}], "]"}], " ", 
      RowBox[{"Cos", "[", 
       RowBox[{
        RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
        RowBox[{"2", " ", "\[Pi]", " ", "y"}]}], "]"}], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"2", " ", "\[Pi]", " ", 
         RowBox[{"Cos", "[", 
          RowBox[{
           RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
           RowBox[{"2", " ", "\[Pi]", " ", "x"}]}], "]"}]}], "-", 
        RowBox[{"4", " ", 
         SuperscriptBox["\[Pi]", "2"], " ", 
         RowBox[{"Cos", "[", 
          RowBox[{"2", " ", "\[Pi]", " ", "x"}], "]"}], " ", 
         RowBox[{"Sin", "[", 
          RowBox[{
           RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
           RowBox[{"2", " ", "\[Pi]", " ", "y"}]}], "]"}]}]}], ")"}]}], "+", 
     RowBox[{"4", " ", 
      SuperscriptBox["\[Pi]", "2"], " ", 
      RowBox[{"Cos", "[", 
       RowBox[{
        RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
        RowBox[{"2", " ", "\[Pi]", " ", "y"}]}], "]"}], " ", 
      RowBox[{"Sin", "[", 
       RowBox[{"2", " ", "\[Pi]", " ", "x"}], "]"}], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", 
         RowBox[{"Sin", "[", 
          RowBox[{
           RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
           RowBox[{"2", " ", "\[Pi]", " ", "x"}]}], "]"}]}], "-", 
        RowBox[{"2", " ", "\[Pi]", " ", 
         RowBox[{"Sin", "[", 
          RowBox[{"2", " ", "\[Pi]", " ", "x"}], "]"}], " ", 
         RowBox[{"Sin", "[", 
          RowBox[{
           RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
           RowBox[{"2", " ", "\[Pi]", " ", "y"}]}], "]"}]}]}], ")"}]}]}], 
    ")"}], " ", 
   RowBox[{"(", 
    RowBox[{"\[Rho]0", "+", 
     RowBox[{
      FractionBox["1", "2"], " ", "\[Rho]1", " ", 
      RowBox[{"(", 
       RowBox[{"1", "+", 
        RowBox[{"Tanh", "[", 
         FractionBox[
          RowBox[{"0.1`", "\[VeryThinSpace]", "-", 
           SqrtBox[
            RowBox[{
             SuperscriptBox[
              RowBox[{"(", 
               RowBox[{
                RowBox[{"-", "0.5`"}], "+", "x"}], ")"}], "2"], "+", 
             SuperscriptBox[
              RowBox[{"(", 
               RowBox[{
                RowBox[{"-", "0.5`"}], "+", "y"}], ")"}], "2"]}]]}], 
          "\[Delta]"], "]"}]}], ")"}]}]}], ")"}]}]}]], "Output",
 CellChangeTimes->{{3.7137151835747633`*^9, 3.713715197281507*^9}, 
   3.713715377068603*^9, {3.713715421668744*^9, 3.713715442688013*^9}, {
   3.713718300209652*^9, 3.713718303371324*^9}, {3.713719094128516*^9, 
   3.7137191371153383`*^9}, 3.7137204530821238`*^9, 3.715696505354093*^9, 
   3.715969644691567*^9, 3.716904709351152*^9, 3.716904766902685*^9, 
   3.717354248157126*^9, {3.71735432615724*^9, 3.717354346852295*^9}, 
   3.717358321406577*^9, 3.71744873648444*^9, 3.717791082128022*^9, 
   3.71779136874577*^9, 3.7177916176324883`*^9, 3.7178065650079308`*^9, 
   3.717857654516314*^9, 3.7294662537271423`*^9, 3.729466675783133*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Periodic", " ", 
    RowBox[{"forces", "?"}]}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"fx", "[", 
     RowBox[{"0", ",", "y", ",", "t"}], "]"}], "-", 
    RowBox[{"fx", "[", 
     RowBox[{"1", ",", "y", ",", "t"}], "]"}]}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"fy", "[", 
     RowBox[{"0", ",", "y", ",", "t"}], "]"}], "-", 
    RowBox[{"fy", "[", 
     RowBox[{"1", ",", "y", ",", "t"}], "]"}]}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"fx", "[", 
     RowBox[{"x", ",", "0", ",", "t"}], "]"}], "-", 
    RowBox[{"fx", "[", 
     RowBox[{"x", ",", "0", ",", "t"}], "]"}]}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"fy", "[", 
     RowBox[{"x", ",", "1", ",", "t"}], "]"}], "-", 
    RowBox[{"fy", "[", 
     RowBox[{"x", ",", "1", ",", "t"}], "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.715696511935528*^9, 3.715696569506382*^9}}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{{3.715696544609459*^9, 3.715696571215577*^9}, 
   3.715969647713558*^9, 3.716904723702054*^9, 3.716904769988472*^9, {
   3.717354336793435*^9, 3.7173543491187983`*^9}, 3.717358326244808*^9, 
   3.717448738550583*^9, 3.717791378260398*^9, 3.717806579395917*^9, 
   3.717857657590087*^9, 3.729466256976911*^9, 3.729466680576256*^9}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{{3.715696544609459*^9, 3.715696571215577*^9}, 
   3.715969647713558*^9, 3.716904723702054*^9, 3.716904769988472*^9, {
   3.717354336793435*^9, 3.7173543491187983`*^9}, 3.717358326244808*^9, 
   3.717448738550583*^9, 3.717791378260398*^9, 3.717806579395917*^9, 
   3.717857657590087*^9, 3.729466256976911*^9, 3.729466680588397*^9}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{{3.715696544609459*^9, 3.715696571215577*^9}, 
   3.715969647713558*^9, 3.716904723702054*^9, 3.716904769988472*^9, {
   3.717354336793435*^9, 3.7173543491187983`*^9}, 3.717358326244808*^9, 
   3.717448738550583*^9, 3.717791378260398*^9, 3.717806579395917*^9, 
   3.717857657590087*^9, 3.729466256976911*^9, 3.729466680596045*^9}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{{3.715696544609459*^9, 3.715696571215577*^9}, 
   3.715969647713558*^9, 3.716904723702054*^9, 3.716904769988472*^9, {
   3.717354336793435*^9, 3.7173543491187983`*^9}, 3.717358326244808*^9, 
   3.717448738550583*^9, 3.717791378260398*^9, 3.717806579395917*^9, 
   3.717857657590087*^9, 3.729466256976911*^9, 3.729466680605584*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"Manipulate", "[", 
  RowBox[{
   RowBox[{"ContourPlot", "[", 
    RowBox[{
     RowBox[{"v", "[", 
      RowBox[{"x", ",", "y", ",", "t"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"x", ",", "0", ",", "1"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"y", ",", "0", ",", "1"}], "}"}], ",", 
     RowBox[{"Contours", "\[Rule]", "10"}]}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"t", ",", "0", ",", "0.001"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.713717913975871*^9, 3.713717944524416*^9}, {
  3.713717998893046*^9, 3.713718046227118*^9}, {3.713718325525816*^9, 
  3.7137183256556396`*^9}, {3.7137183565887938`*^9, 3.7137183566716137`*^9}, {
  3.713718842734448*^9, 3.7137188710275707`*^9}, {3.7137237118451653`*^9, 
  3.713723714042118*^9}, {3.729466284418435*^9, 3.7294662971641397`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"\[Rho]0", " ", "=", " ", "1.0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Rho]1", " ", "=", " ", 
   RowBox[{
    RowBox[{"-", "1.0"}], "+", "1000"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Delta]", " ", "=", " ", "0.01"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"ContourPlot", "[", 
  RowBox[{
   RowBox[{"\[Rho]", "[", 
    RowBox[{"x", ",", "y", ",", "0"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", "0", ",", "1"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"y", ",", "0", ",", "1"}], "}"}], ",", 
   RowBox[{"Contours", "->", "20"}]}], "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Rho]", "[", 
   RowBox[{"0.5", ",", "0.5", ",", "0"}], "]"}], "/", 
  RowBox[{"\[Rho]", "[", 
   RowBox[{"1", ",", "1", ",", "0"}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.71596965882487*^9, 3.715969785676296*^9}, {
  3.7159698210126047`*^9, 3.715969854032394*^9}, {3.715970539524322*^9, 
  3.715970575472484*^9}, {3.715970735131856*^9, 3.7159707419375134`*^9}, {
  3.715978463957518*^9, 3.7159784651295137`*^9}}],

Cell[BoxData["999.9999979409074`"], "Output",
 CellChangeTimes->{
  3.715969671681885*^9, {3.7159697023931923`*^9, 3.71596972850812*^9}, {
   3.715969770130165*^9, 3.715969774416934*^9}, {3.715969823216819*^9, 
   3.715969868117879*^9}, {3.715970541241802*^9, 3.715970575920765*^9}, 
   3.715970742476511*^9, 3.715978466930719*^9, 3.7294666856438*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"D", "[", 
   RowBox[{
    RowBox[{"p", "[", 
     RowBox[{"x", ",", "y", ",", "t"}], "]"}], ",", "x"}], "]"}], "//", 
  "FullSimplify"}]], "Input",
 CellChangeTimes->{{3.716904802372607*^9, 3.716904803365445*^9}}],

Cell[BoxData[
 RowBox[{"4", " ", 
  SuperscriptBox["\[Pi]", "2"], " ", 
  RowBox[{"Cos", "[", 
   RowBox[{"2", " ", "\[Pi]", " ", 
    RowBox[{"(", 
     RowBox[{"t", "-", "x"}], ")"}]}], "]"}], " ", 
  RowBox[{"Cos", "[", 
   RowBox[{"2", " ", "\[Pi]", " ", 
    RowBox[{"(", 
     RowBox[{"t", "-", "y"}], ")"}]}], "]"}]}]], "Output",
 CellChangeTimes->{{3.716904798182521*^9, 3.716904803801478*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"D", "[", 
   RowBox[{
    RowBox[{"p", "[", 
     RowBox[{"x", ",", "y", ",", "t"}], "]"}], ",", "y"}], "]"}], "//", 
  "FullSimplify"}]], "Input",
 CellChangeTimes->{{3.7169048001843233`*^9, 3.7169048109077377`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", "4"}], " ", 
  SuperscriptBox["\[Pi]", "2"], " ", 
  RowBox[{"Sin", "[", 
   RowBox[{"2", " ", "\[Pi]", " ", 
    RowBox[{"(", 
     RowBox[{"t", "-", "x"}], ")"}]}], "]"}], " ", 
  RowBox[{"Sin", "[", 
   RowBox[{"2", " ", "\[Pi]", " ", 
    RowBox[{"(", 
     RowBox[{"t", "-", "y"}], ")"}]}], "]"}]}]], "Output",
 CellChangeTimes->{{3.716904800466488*^9, 3.716904811201728*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"CellPrint", "[", 
  RowBox[{"[", 
   RowBox[{
    RowBox[{"HoldForm", "[", 
     RowBox[{"fx", "[", 
      RowBox[{"x", ",", "y", ",", "t"}], "]"}], "]"}], ",", "\"\<Output\>\"", 
    ",", 
    RowBox[{"AutoMultiplicationSymbol", "\[Rule]", "True"}]}], "]"}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.7173548332167797`*^9, 3.717354904647806*^9}}],

Cell[BoxData[
 RowBox[{"CellPrint", "\[LeftDoubleBracket]", 
  RowBox[{
   TagBox[
    RowBox[{"fx", "[", 
     RowBox[{"x", ",", "y", ",", "t"}], "]"}],
    HoldForm], ",", "\<\"Output\"\>", ",", 
   RowBox[{"AutoMultiplicationSymbol", "\[Rule]", "True"}]}], 
  "\[RightDoubleBracket]"}]], "Output",
 CellChangeTimes->{3.717354852169023*^9, 3.7173549152477007`*^9},
 AutoMultiplicationSymbol->True]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"fy", "[", 
   RowBox[{"x", ",", "y", ",", "t"}], "]"}], "//", "InputForm"}]], "Input",
 CellChangeTimes->{{3.717355013352757*^9, 3.717355014922484*^9}, {
  3.717355057384183*^9, 3.71735505745322*^9}, {3.717791745635724*^9, 
  3.717791753211072*^9}, {3.7178066786883087`*^9, 3.7178066789195833`*^9}, {
  3.717806716945622*^9, 3.7178067179080467`*^9}}],

Cell["\<\
-4*Pi^2*Sin[2*Pi*t - 2*Pi*x]*Sin[2*Pi*t - 2*Pi*y] - 
 2*(2*Pi^2*\[Mu]1*Cos[2*Pi*x]*Cos[2*Pi*t - 2*Pi*x]*Cos[2*Pi*y] + 
   2*Pi^2*(\[Mu]0 + \[Mu]1 + \[Mu]1*Cos[2*Pi*y]*Sin[2*Pi*x])*Sin[2*Pi*t - \
2*Pi*x] - 
   8*Pi^3*\[Mu]1*Cos[2*Pi*t - 2*Pi*y]*Sin[2*Pi*x]^2*Sin[2*Pi*y] + 
   8*Pi^3*Sin[2*Pi*x]*(\[Mu]0 + \[Mu]1 + \
\[Mu]1*Cos[2*Pi*y]*Sin[2*Pi*x])*Sin[2*Pi*t - 2*Pi*y]) + 
 (-2*Pi*Cos[2*Pi*t - 2*Pi*x] - 4*Pi^2*Cos[2*Pi*t - 2*Pi*y]*Sin[2*Pi*x] + 
   2*Pi*Cos[2*Pi*x]*Cos[2*Pi*t - 2*Pi*y]*(2*Pi*Cos[2*Pi*t - 2*Pi*x] - 
     4*Pi^2*Cos[2*Pi*x]*Sin[2*Pi*t - 2*Pi*y]) + 4*Pi^2*Cos[2*Pi*t - \
2*Pi*y]*Sin[2*Pi*x]*
    (-Sin[2*Pi*t - 2*Pi*x] - 2*Pi*Sin[2*Pi*x]*Sin[2*Pi*t - 2*Pi*y]))*
  (\[Rho]0 + (\[Rho]1*(1 + Tanh[(0.1 - Sqrt[(-0.5 + x)^2 + (-0.5 + y)^2])/\
\[Delta]]))/2)\
\>", "Output",
 CellChangeTimes->{
  3.717355015260579*^9, 3.717355057957321*^9, {3.717791746216078*^9, 
   3.717791753461899*^9}, 3.717806679157455*^9, 3.717806718151379*^9, 
   3.717806756830654*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"\[Rho]0", " ", "=", " ", "1.0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Rho]1", " ", "=", " ", 
   RowBox[{"1.0", "+", "10"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Mu]0", " ", "=", " ", "0.001"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Mu]1", " ", "=", " ", 
   RowBox[{
    RowBox[{"-", "\[Mu]0"}], "+", "1"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Delta]", " ", "=", " ", "0.01"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Manipulate", "[", 
  RowBox[{
   RowBox[{"ContourPlot", "[", 
    RowBox[{
     RowBox[{"fy", "[", 
      RowBox[{"x", ",", "y", ",", "t"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"x", ",", "0", ",", "1"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"y", ",", "0", ",", "1"}], "}"}], ",", 
     RowBox[{"Contours", "\[Rule]", "10"}], ",", 
     RowBox[{"PlotLegends", "\[Rule]", "Automatic"}]}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"t", ",", "0", ",", "0.1"}], "}"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.717358156058093*^9, 3.717358245129712*^9}, {
  3.7174487594454203`*^9, 3.717448784062738*^9}, {3.71744887398997*^9, 
  3.7174488774384413`*^9}, {3.729466700372604*^9, 3.729466700728937*^9}, {
  3.729467218057976*^9, 3.72946723459153*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "Tangential", " ", "traction", " ", "at", " ", "top", " ", "and", " ", 
    "bottom", " ", "face"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"\[Mu]", "[", 
    RowBox[{"x", ",", "y", ",", "t"}], "]"}], 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"D", "[", 
      RowBox[{
       RowBox[{"u", "[", 
        RowBox[{"x", ",", "y", ",", "t"}], "]"}], ",", "y"}], "]"}], "+", " ", 
     RowBox[{"D", "[", 
      RowBox[{
       RowBox[{"v", "[", 
        RowBox[{"x", ",", "y", ",", "t"}], "]"}], ",", "x"}], "]"}]}], 
    ")"}]}]}]], "Input",
 CellChangeTimes->{{3.7177910183756227`*^9, 3.717791074095456*^9}, {
  3.71779125004282*^9, 3.7177912581423492`*^9}}],

Cell[BoxData[
 RowBox[{"2", " ", "\[Pi]", " ", 
  RowBox[{"Cos", "[", 
   RowBox[{
    RowBox[{"2", " ", "\[Pi]", " ", "t"}], "-", 
    RowBox[{"2", " ", "\[Pi]", " ", "x"}]}], "]"}], " ", 
  RowBox[{"(", 
   RowBox[{"\[Mu]0", "+", "\[Mu]1", "+", 
    RowBox[{"\[Mu]1", " ", 
     RowBox[{"Cos", "[", 
      RowBox[{"2", " ", "\[Pi]", " ", "y"}], "]"}], " ", 
     RowBox[{"Sin", "[", 
      RowBox[{"2", " ", "\[Pi]", " ", "x"}], "]"}]}]}], ")"}]}]], "Output",
 CellChangeTimes->{{3.717791074728895*^9, 3.717791088323224*^9}, {
   3.717791246876535*^9, 3.7177912585491247`*^9}, {3.717791336071959*^9, 
   3.717791350906612*^9}, 3.717791385245235*^9, 3.717792645806684*^9, 
   3.717806589758082*^9, 3.717857664600211*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "Normal", " ", "traction", " ", "at", " ", "top", " ", "and", " ", 
    "bottom", " ", "face"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"-", 
     RowBox[{"p", "[", 
      RowBox[{"x", ",", "y", ",", "t"}], "]"}]}], "+", 
    RowBox[{"2", 
     RowBox[{"\[Mu]", "[", 
      RowBox[{"x", ",", "y", ",", "t"}], "]"}], 
     RowBox[{"D", "[", 
      RowBox[{
       RowBox[{"v", "[", 
        RowBox[{"x", ",", "y", ",", "t"}], "]"}], ",", "y"}], "]"}]}]}], "//",
    "InputForm"}]}]], "Input",
 CellChangeTimes->{{3.71780899757131*^9, 3.7178090216842327`*^9}, {
  3.717809098251532*^9, 3.7178091076245127`*^9}, {3.717809201774146*^9, 
  3.717809228762733*^9}}],

Cell["\<\
8*Pi^2*Cos[2*Pi*t - 2*Pi*y]*Sin[2*Pi*x]*(\[Mu]0 + \[Mu]1 + \
\[Mu]1*Cos[2*Pi*y]*Sin[2*Pi*x]) + 
 2*Pi*Cos[2*Pi*t - 2*Pi*y]*Sin[2*Pi*t - 2*Pi*x]\
\>", "Output",
 CellChangeTimes->{{3.717809104244463*^9, 3.717809109406019*^9}, {
   3.717809202593155*^9, 3.717809229111446*^9}, 3.7178576660192547`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"trac", "[", 
   RowBox[{"x_", ",", "y_", ",", "t_"}], "]"}], "=", 
  RowBox[{
   RowBox[{"-", 
    RowBox[{"p", "[", 
     RowBox[{"x", ",", "y", ",", "t"}], "]"}]}], "+", 
   RowBox[{"2", 
    RowBox[{"\[Mu]", "[", 
     RowBox[{"x", ",", "y", ",", "t"}], "]"}], 
    RowBox[{"D", "[", 
     RowBox[{
      RowBox[{"v", "[", 
       RowBox[{"x", ",", "y", ",", "t"}], "]"}], ",", "y"}], 
     "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.717852484006867*^9, 3.717852489067326*^9}, {
  3.717852562818201*^9, 3.717852565223977*^9}, {3.7178576292823553`*^9, 
  3.717857630359495*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", 
   RowBox[{"p", "[", 
    RowBox[{"x", ",", "y", ",", "t"}], "]"}]}], "+", 
  RowBox[{"2", " ", 
   RowBox[{"\[Mu]", "[", 
    RowBox[{"x", ",", "y", ",", "t"}], "]"}], " ", 
   RowBox[{
    SuperscriptBox["v", 
     TagBox[
      RowBox[{"(", 
       RowBox[{"0", ",", "1", ",", "0"}], ")"}],
      Derivative],
     MultilineFunction->None], "[", 
    RowBox[{"x", ",", "y", ",", "t"}], "]"}]}]}]], "Output",
 CellChangeTimes->{
  3.717852489804448*^9, 3.717852529671657*^9, 3.717852567525247*^9, {
   3.717857630819322*^9, 3.7178576445569477`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"trac", "[", 
   RowBox[{"0", ",", "1", ",", "t"}], "]"}], "-", 
  RowBox[{"trac", "[", 
   RowBox[{"1", ",", "1", ",", "t"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.7178524905294247`*^9, 3.7178525892724247`*^9}}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{{3.717852578614255*^9, 3.717852589621361*^9}, 
   3.7178576321350603`*^9}]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.7178525836059*^9, 3.717852584903145*^9}}]
},
WindowSize->{800, 675},
WindowMargins->{{83, Automatic}, {61, Automatic}},
FrontEndVersion->"10.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (June 27, \
2014)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 76, 1, 28, "Input"],
Cell[637, 23, 3934, 102, 182, "Input"],
Cell[4574, 127, 1187, 38, 63, "Input"],
Cell[CellGroupData[{
Cell[5786, 169, 644, 17, 46, "Input"],
Cell[6433, 188, 620, 9, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[7090, 202, 4119, 119, 131, "Input"],
Cell[11212, 323, 4964, 127, 179, "Output"],
Cell[16179, 452, 6071, 155, 221, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[22287, 612, 943, 26, 97, "Input"],
Cell[23233, 640, 378, 5, 28, "Output"],
Cell[23614, 647, 378, 5, 28, "Output"],
Cell[23995, 654, 378, 5, 28, "Output"],
Cell[24376, 661, 378, 5, 28, "Output"]
}, Open  ]],
Cell[24769, 669, 833, 18, 46, "Input"],
Cell[CellGroupData[{
Cell[25627, 691, 1088, 26, 97, "Input"],
Cell[26718, 719, 351, 5, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[27106, 729, 248, 7, 28, "Input"],
Cell[27357, 738, 402, 11, 32, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[27796, 754, 252, 7, 28, "Input"],
Cell[28051, 763, 420, 12, 32, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[28508, 780, 371, 10, 28, "Input"],
Cell[28882, 792, 399, 10, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[29318, 807, 385, 7, 28, "Input"],
Cell[29706, 816, 982, 19, 164, "Output"]
}, Open  ]],
Cell[30703, 838, 1293, 31, 131, "Input"],
Cell[CellGroupData[{
Cell[32021, 873, 727, 21, 46, "Input"],
Cell[32751, 896, 721, 16, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[33509, 917, 742, 21, 46, "Input"],
Cell[34254, 940, 309, 6, 44, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[34600, 951, 617, 18, 28, "Input"],
Cell[35220, 971, 591, 18, 32, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[35848, 994, 249, 6, 28, "Input"],
Cell[36100, 1002, 122, 2, 28, "Output"]
}, Open  ]],
Cell[36237, 1007, 90, 1, 28, "Input"]
}
]
*)

(* End of internal cache information *)
